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Abstract 

Using the recent measurement of the total production rate for Bg and -B* mesons 
in electron-positron annihilation to determine the strange quark mass parameter in the 
b ^ Bs, B* fragmentation functions we calculate the momentum distributions of the 
-Bo and Bt mesons. 
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Recently the probability for a heavy quark to hadronize into a strange B or D meson 
has been measured at LEP by the DELPHI Collaboration These data can be used to 
determine the strange quark mass parameter, m^, in the heavy quark fragmentation functions 
and, therefore, to predict the momentum distributions of the Bg and B* mesons produced 
in the fragmentation process. Using the heavy quark fragmentation functions calculated 
in perturbative QCD at the scale of the heavy quark for the initial condition we have 
numerically integrated the Altarelli-Parisi equation to obtain the momentum distributions 
of the Bg and B* mesons at the scale of the Z mass. The resulting distributions are only 
moderately sensitive to the value of m^, which we find to be about 300 MeV/c^. 

The fragmentation function for the process b ^ Bg at the scale of the b mass is easily 
obtained from the b ^ Be fragmentation functions calculated using perturbative QCD by 
Braaten, Cheung, and Yuan [Q: 
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Here z is the fraction of the b quark momentum carried by the Bg meson, r = ms/{mb + mg), 
and -R(O) is the Bg meson S-wave radial wavefunction at the origin. The corresponding 
fragmentation function for b ^ B* is 
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The wavefunction at the origin, -R(O), can be determined from the Bg meson decay constant, 
fs^, from the Van Royen-Weisskopf [Q relation, modified for color: 
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In Eqs. (|T|) and (|^) the running coupling constant ^^(/i) = as{Mz) / {l+boas{Mz) log(/i/M^)), 
where 60 = (33 — 2nf)/6iT, Uf is the number of flavors at the scale /i, as{Mz) = 0.12, and 
we have chosen the scale to be 2ms • The scale /iq of the fragmentation functions is of the 
order of and we choose it to be + 2ms [0- Integrating the Altarelli-Parisi equation 
for the fragmentation functions from the initial scale fiQ to the Z mass gives the momentum 
distributions of the Bg and B* mesons produced by the fragmentation process in e+e^ anni- 
hilation at the Z resonance. To determine the strange quark mass parameter, m^, we use the 
recent measurement by the DELPHI collaboration [|l| of the probability, f^, that a weakly 
decaying strange heavy meson is produced during the hadronization of a heavy quark. In 
addition to the direct production b ^ Bs, the Bs meson can also result from the decays of 
excited states of Bs and Be mesons. However, these have been excluded in the measurement 
of f^, except for the radiative decay of the first excitation: B* ^ Bg + 7. Since the Bs — B* 
mass difference is so small we can determine from the integrated fragmentation functions, 
which gives the measured probability f^: 

f: = f' dz [Di^b. {z, fio) + Di^Bt /"o)] . (4) 
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Since the integrated fragmentation functions are independent of the scale, we can choose 
the scale to be /io in Eq. (Q) and first determine and then calculate the distributions 
Di^B^{z, Mz) and Di_^b*{z, Mz) at the scale of the Z mass by numerically integrating the 
Altarelli-Parisi equation. As input we fix the heavy b quark mass at m^ = 5.0 GeV/c^. To 
determine i?(0) we use Eq. (g) with the value fs, = 207 ± 34 ± 22 MeV from the lattice 
calculations of Bernard, Labrenz, and Soni 0. Then from the measured value 

= 0.19 ± 0.06 ± 0.08 (5) 

we find (adding the two errors in quadrature) 

ms = 318 _24 MeV , (6) 

which is in a reasonable range 0. The estimated error in reflects only the experimental 
uncertainty in f^. 



Using this value of rris we have calculated the momentum distributions Di^^^{z) and 
Di^Qt{z), which are shown in Figs. |l|. Figure |l](a) shows the fragmentation functions at the 
initial scale /io = mi, + 2ms, while Fig. ^(b) shows the fragmentation functions at the scale 
/i = Mz- We have verified that the shapes of these distributions are only moderately sensitive 
to the value of m^, although the peak does shift towards larger values of z if is decreased, 
as one might expect. In fact, the maxima in both Di^^^^z, Mz) and Di_^q*{z, Mz) occur at 

B:) = 0.92 . (7) 

The average momentum fraction (z) at the Z mass scale is also not too different for the Bg 
and B* contributions to the spectrum: 

= 0.65 and (z)-,^^. = 0.68 . (8) 

The spectrum also shows some evidence for approximate heavy quark spin symmetry at 
the level of about 50%. The ratio / dzDi_^^,/ J dzDi^^^ is about 2.1 compared to 3, as 
expected in the heavy quark limit from the spin independence. In terms of the ratio Py = 
B*/{Bs + B*), which heavy quark symmetry predicts to be 0.75, we find 0.68. To compare 
the b ^ Bg and b ^ B* spectra individually with experiments, the two contributions b ^ Bg 
and b ^ B* ^ Bg^ have to be experimentally separated by observing the rather low energy 
photon, which has so far not been done. 

To explore the sensitivity of our calculation to the assumptions we also considered the 
fragmentation of a charm quark into a Dg or D* meson. Following ref. we repeated the 
calculation assuming the probability that a Dg or D* meson is produced in the hadronization 
of a c quark is the same as the probability that a Bg or B* meson is produced in the 
hadronization of a 6 quark; that is, is independent of the heavy quark flavor. Using the 
value = 0.19 ± 0.06 ± 0.08 0, the lattice calculation Jd, = 230 ± 30 ± 18 MeV |], and 
choosing rric = 1.5 GeV/c^ 0, we found 

m, = 267^?7 MeV, (9) 
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which is in reasonable agreement with the value of obtained above in the case. With 
this value of we also calculated the momentum spectra for c Dg, for c ^ D*, and their 
sum at the initial scale /Xq = + 2ms as well as at the scale fi = Mz- They are shown, 
respectively, in Figs, ^(a) and 0(b). As expected, the Ds spectra are softer than the Bg 
spectra shown in Figs. 0(a) and ^(b) because of the smaller charm quark mass. The peaks 
of the spectra occur at 0.68 for Dg and at 0.75 for D* and the mean momentum fractions 
are {z) = 0.49 for Dg and (z) = 0.52 for L>* at the scale n = Mz [Fig. |(b)]. The ratio of 
the two fragmentation probabilities for the Dg case also differs significantly more from the 
heavy quark spin symmetry prediction of 3, being about 1.6, than in the Bg case, as one 
expects. Also, the ratio Py = D*g/{Ds + D*g) is 0.62, rather than 0.75. 

To further test the sensitivity of our results to the input parameters assumed we have 
repeated the calculation of the Dg and D*g spectra using the measurement of the decay 
constant fo^ recently reported by the CLEO collaboration 0. Although this measured value 
/d, = 344±37±52±42 MeV differs substantially from the lattice result /d, = 230±30±18 
MeV 0, we found 

mg = 306 Itl MeV , (10) 

which is not significantly different from the value above [Eq. (|^)], and the changes in the 
fragmentation spectra shown in Fig. ^ are negligible. 

To summarize, we have calculated the momentum distributions of Bg and B* mesons 
produced by b quark fragmentation. The measured total production probability was used 
to determine the strange quark mass parameter. The result for nig was quite reasonable, 
being well above the current quark mass value, but lower than values often used for the 
constituent quark mass. We note that at the scale 2mg that as{2'mg) ^ 0.8 is uncomfortably 
large, which could mean there are significant corrections, which have not been calculated, to 
the fragmentation function at the heavy quark scale, which was calculated using perturbative 
QCD. But this only affects the normalization and we determined the parameter rrig from the 
experimental value of the total production probability. Indeed, we repeated the calculation 
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for Dg and D* meson production by fragmentation and found substantially the same value 
for nis and verified that the momentum distributions were only moderately sensitive to the 
value of rUs- The results were also not sensitive to the values assumed for the heavy quark 
masses or the meson decay constants. Essentially, the experimental production probability 
fixes the normalization, the perturbative QCD result determines the shape of the momen- 
tum spectrum at the heavy quark mass scale, and the Altarelli-Parisi equation governs the 
evolution up to the scale Mz] so perhaps it is not so surprising that the results are not too 
sensitive to the particular values of the parameters. 

This work was supported by the U. S. Department of Energy, Division of High Energy 
Physics, under Grant DE-FG02-91-ER40684. 
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The first two errors are the statistical and systematic errors, respectively, while the last 
error reflects the uncertainty in the absolute ^tt"*" branching ratio. 
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Figure Captions 

1. (a) Momentum distributions for the fragmentation processes b ^ Bg, b ^ B*, and 

their sum at the scale = n^b + 2ms. 

(b) Momentum distributions for the fragmentation processes b ^ Bs, b ^ B*, and 
their sum at the scale = Mz- 

2. (a) Momentum distributions for the fragmentation processes c — > Dg, c D*, and 
their sum at the scale /Xq = rric + 2ms. 

(b) Momentum distributions for the fragmentation processes c Dg, c —>■ D*, and 
their sum at the scale /i — Mz. 
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This figure "figl-l.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9406207vl 



This figure "figl-2.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9406207vl 



